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We present new constraints on cosmic variations of Newton's gravitational constant by making 
use of the latest CMB data from WMAP, BOOMERANG, CBI and ACBAR experiments and 
independent constraints coming from Big Bang Nucleosynthesis. We found that current CMB data 
provide constraints at the ~ 10% level, that can be improved to ~ 3% by including BBN data. We 
show that future data expected from the Planck satellite could constrain G at the ~ 1.5% level 
while an ultimate, cosmic variance limited, CMB experiment could reach a precision of about 0.4%, 
competitive with current laboratory measurements. 
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I. INTRODUCTION 

Since Cavendish's first measurement in 1798 (jjj), 
Newton's Gravitational constant remains one of the most 
elusive constants in physics. The past two decades did 
not succeed in substantially improving our knowledge of 
its value from the precision of 0.05% reached in 1942 
(see H). To the contrary, the variation between different 
measurements forced the CODATA committee 1 , which 
determines the internationally accepted standard values, 
to increase the uncertainty from 0.013% for the value 
quoted in 1987 to the one order of magnitude larger un- 
certainty of 0.15% for the 1998 "official" value (g). Re- 
cent laboratory measurements (see e.g. g) point towards 
an uncertainty at the level of ~ 0.4%, while other works 
claim an improved precisions below 0.01% (||). Analysis 
of the secular variation of the period of nonradial pulsa- 
tions of the white dwarf G117-B15A (0) has produced 
complementary constraints at ~ 0.1% level. 

Measurements of the Cosmic Microwave Background 
(CMB, hereafter) temperature and polarization 
anisotropy have been suggested as a possible tool 
for determining the value of G (see @). In recent years, 
CMB temperature and polarization anisotropy have 
been measured with great precision from experiments 
as WMAP (fH), BOOMERANG (JkJ), CBI gl| and 
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ACBAR (JL2[). The impressive agreement between those 
measurements and the expectations of the standard 
model of structure formation have paved the way to 
the use of cosmology as a new laboratory where to 
test physical hypothesis at energies and scales not 
reachable on earth. Since a variation in G affects CMB 
temperature and polarization anisotropy, changing the 
position and the amplitude of the acoustic peaks present 
in the corresponding angular power spectra, it is indeed 
possible to infer new and independent constraints on G 
from CMB data. 

In this paper we follow this timely line of investigation. 
Respect to previous works (most notably Q) we update 
the CMB constraints on G by using the most recent CMB 
data (most notably, WMAP) and by also including com- 
plementary information from Big Bang Nucleosynthesis 
(hereafter, BBN, see jp| for a complete review). As al- 
ready shown in several papers (see e.g. E3], fL5|), any 
variation in G changes the Hubble parameter at BBN 
given by H ~ y/Gg^T 2 where g* counts the number of 
relativistic particles species and T is the temperature of 
the Universe. Since the predicted amount of light ele- 
ments depends crucially on the comparison between the 
expansion rate H and, for example, the neutron-proton 



conversion rate T„ 



G 2 F T 5 , 



where Gf is the Fermi 
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constant, any change in G can be strongly constrained 
by combining BBN predictions with observations of pri- 
mordial elements. Moreover, we also discuss the ability of 
next CMB experiments as Planck ([l6|]) to constrain G, 
including the possibility of a "cosmic variance limited" 
survey. 

Any cosmological constraint is, however, indirect and, 
in the case of the CMB data, depends on the assumed 
theory of structure formation. The major caveat in our 
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case is the assumption of a cosmological constant, or dark 
energy component, the nature of which is puzzling and 
unknown (for a recent review, see e.g. |l7| , [Q, [^9|). 
While the derived constraints will therefore be model de- 
pendent, it is interesting that a major alternative to a 
dark energy component, i.e. modified gravity theories, 
could be parameterised by introducing an effective value 
of Newton's constant G c g , that could not only be differ- 
ent from the local value of G but also spatial and time 
dependent (see e.g. |2(j, |2l),[|2|). Moreover, if dark en- 
ergy interacts with dark matter, there is a change in the 
background evolution of the universe leading to an effec- 
tive G e g for the matter component (see e.g. and to 
a possible change in the cosmic bound on G. 

In this respect, the search for variations in Newton's 
constant using cosmological data could also play a role 
in the understanding of the dark sector. If the New- 
ton's constant inferred from cosmology will turn out to 
be different from the local value, then this may suggest a 
modification of gravity at large scale or a more complex 
interacting dark energy scenario. Since an interacting 
dark energy or a modified gravity theory could be re- 
sponsible for a variation of G in the late universe, we 
also consider the possibility of a rcdshift dependence of 
G. 

Our paper is therefore organized as follows: in the next 
section we briefly describe the effects of a variation in G 
on CMB temperature and polarization anisotropy. In 
Section III we describe our method of analysis and the 
datasets considered. In Section IV we present our results 
and, finally, in Section V we derive our conclusions. 

II. THE IMPACT OF G ON RECOMBINATION 
AND THE CMB 

Following Q] we parameterize the deviations from 
Newton's gravitational constant by introducing a dimcn- 
sionless parameter Xq such that 

G^X G G (1) 

As showed in 0, expressing the perturbed quantities 
in Fourier space, a variation in Newton's gravitational 
constant is equivalent in a simple re-scaling of the wave 
numbers. No preferred cosmological scale is introduced 
by varying G and the density fluctuations produced by 
a mode of wavevector k in a universe with Xq ^ 1 have 
equivalent dynamics of a mode with k' = k/Xc in a uni- 
verse with Xq = f • 

However the physics of recombination does introduce 
a preferred timescale and it will actually change when 
varying Xg ■ This is clearly shown in Figure |l| where the 
ionization fraction x e at different redshift z, computed 
with a modified version of RECFAST ||, is plotted for 
different values of Xq- The ionization fraction x e is just 
the free electron number density n e divided by the total 
number density of hydrogen nuclei (free and bound) uh ■ 
As we can see, higher (lower) values of the gravitational 



constant yields a delayed (accelerated) period of recombi- 
nation. A change in the number density of free electrons 
n e in function of the conformal time r, changes the vis- 
ibility function g(r), written in terms of the opacity for 
Thomson scattering k as 

g(r) = /iexp(— n) = — d/dr exp(— k) (2) 

with 

/TO 
an e (r)dr, (3) 

where ut is the Thomson scattering cross section, a is 
the scale factor and k = GTan e . 

This clearly affects the CMB temperature anisotropy 
that can be written as an integral along the line of sight 
over sources, 

AT(n, k) = f ° dr S(k, r)e 4k ' fiI) ( r ).g(r) (4) 
Jo 

where S(k,r) is the anisotropy source term (see [Q) and 
D(t) is the distance from the observer to a point along 
the line of sight at conformal time r. 
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FIG. 1: Ionization fraction in function of redshift for different 
values of Xg- 

In Figure |^ we plot the CMB temperature and polar- 
ization spectra computed from a modified version of the 
CAMB j30| code. The effect of modified recombination 
is clear. Namely, varying Xg changes the recombination 
process, shifting g(r) along the conformal time r. The 
net effect is a damping or enhancement of the acoustic os- 
cillations and a shift of the Doppler peaks in the angular 
scales. This mechanism could mimic an extra injection 
or absorption of Lyman-a photons at last scattering, as 
already analyzed in several recent papers (see e.g. p5|), 
and it would be difficult to disentangle the two scenarios. 

Another important aspect to consider is a possible red- 
shift dependence of G. If interacting dark energy or a 
modification to general relativity are responsible for the 
current accelerated expansion of the universe, it is indeed 
possible that this could result in an observed cosmic value 
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FIG. 2: From Top to Bottom: Temperature, Polarization and 
cross Temperature-Polarization power spectra in function of 
variations in Ag- 

of G different from the one obtained from local measure- 
ments. Moreover, it is plausible to think that this kind 
of deviation of G will be triggered by acceleration, i.e., 
to be conservative, will appear at redshift 0.1 < z < 2. 

We have therefore considered two possible parameter- 
izations for a redshift-depcndcnt gravitational constant. 
A first parameterization, that somewhat ties the change 
in G with the appearance of dark energy is to consider: 

G(z) = G + AG(1 — a) (5) 

where the variation AG is equal to G(Xq — 1). This pa- 
rameterization, similar to the one proposed in p7| for 
the dark energy equation of state, has the advantage of a 



smooth transition between the value of G today to XqG 
in the past, when z » 1. However the redshift of tran- 
sition between these two values is not an independent 
variable. 

We have therefore considered a second possible param- 
eterization as: 

G(z) = G[l-(l-\ 2 G )H(z-z t )} (6) 

where H{x) is the Heaviside function (H(x) for x < and 
H (x) = 1 for x > 0) and is the redshift of transition 
between the two values (local and past) of G. 
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FIG. 3: Temperature power spectrum when the gravitational 
constant varies such that Ag = 0.9. The graph shows the 
effects of the smooth transition parameterization described in 
equation |H| (green line) and of the Heaviside parameterization 
of equation |fj for different redshifts of transition zt between 
0.1 and 2. 

In Figure |^ we plot different power spectra computed 
considering the two parameterizations described using a 
fixed value of Xq = 0.9. As we can see, introducing 
a redshift dependent variation in G increases the CMB 
anisotropy at large angular scales. On sub-Hubble scales, 
the Einstein equations in an expanding space-time reduce 
to the Poisson equation 

A$ = AitGpa 2 8 (7) 

that relates the gravitational potential <I> to the density 
contrast 5. If a redshift variation in G occurs, this will 
clearly change the gravitational potential, the density 
growth function and large scale CMB anisotropy through 
the Integrated Sachs Wolfe effect (ISW hereafter, see e.g. 
p6|). Since a large ISW signal is at odds with current 
WMAP data, a varying with redshift G is strongly con- 
strained, as we will see in the next section. 

III. ANALYSIS METHOD 

We constrain variations in the Newton's constant with 
current CMB data by making use of the publicly avail- 
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able Markov Chain Monte Carlo package cosmomc p8| . 
Other than Ag we sample the following set of cosmologi- 
cal parameters, adopting flat priors on them: the physical 
baryon and CDM densities, w& = f2f,/i 2 and lo c = f2 c /i 2 , 
the Hubble parameter, Hq, the scalar spectral index, 
n s , the normalization, lnl0 10 A s (fc = 0.05/ Mpc) and the 
reionization optical depth t. 

As discussed in the previous section, we will also con- 
sider the possibility of a variation with rcdshift in G and 
we will consider as extra parameter the redshift of tran- 
sition zt- 

The MCMC convergence diagnostic tests are per- 
formed on 4 chains using the Gelman and Rubin "vari- 
ance of chain mean" /"mean of chain variances" R — 1 
statistic for each parameter. Our 1 — D and 2 — D 
constraints are obtained after marginalization over the 
remaining "nuisance" parameters, again using the pro- 
grams included in the cosmomc package. We use a cos- 
mic age top-hat prior as 10 Gyr < t < 20 Gyr. We 
include the five-year WMAP data |J (temperature and 
polarization) with the routine for computing the likeli- 
hood supplied by the WMAP team (we will refer to this 
analysis as WMAP5). 

Moreover, in order to test the effect of current polariza- 
tion measurements on constraining Ag we also considered 
the combination of the WMAP data with the polariza- 
tion results coming from the BOOMERANG ([gof) and 
CBI (jll|]) experiments. We will refer to this analysis as 
WMAP5+POL. 

Together with the WMAP data we also consider the 
small-scale CMB measurements of ACBAR [Q (wc will 
refer to this analysis as WMAP5+ ACBAR). 

Finally, we forecast future constraints on Ag simulat- 
ing a set of mock data with a fiducial model given by the 
best fit WMAP5 model with Ag = 1 and experimental 
noise described by: 
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where w x l 2 is the temperature noise level (we consider a 
factor V2 larger for polarization noise) and 9 is the beam 
size. We considered two future datasets. The first, based 
on the experimental specifications of the PLANCK SUR- 
VEYOR mission, with w 1/2 = 58uK and 6> FWHM = 7.1' 
equivalent to the 143 GHz channel (see [|l(| ) . The second 
dataset is a cosmic variance limited experiment (CVL 
hereafter) with no experimental noise for both tempera- 
ture and polarization anisotropy and £ max = 2500. 

Constraints on Ag are also computed using standard 
BBN theoretical predictions as provided by the new nu- 
merical code described in [^|(3^], which includes a full 
updating of all rates entering the nuclear chain based 
on the most recent experimental results on nuclear cross 
sections. The BBN predictions are compared with the 
D/H abundance ratio of |3l[] obtained including a new 
measurement in a metal poor damped Lyman-a system 



along the line of sight of QSO SDSS1558-0031 

D/H = (2.82t°J£) ■ 10" 5 (9) 

We use the uncertainty as quoted in |3l| , computed by a 
jackknifc analysis. 

IV. RESULTS 
A. Constant G with redshift 

We report in Table | the constraints obtained on Ag 
analyzing the datasets mentioned in the previous sec- 
tion. As we can see, current CMB data only provide 
a constraints at about ~ 15% level. The WMAP con- 
straint is improved by ~ 10% when temperature and 
polarization anisotropy data from BOOMERANG and 
CBI is included and by ~ 30% when the small scale tem- 
perature angular spectrum data from ACBAR is added. 
However, as we can see from the Table, the major im- 
provement comes from BBN: in this case the constraint 
WMAP+BBN reaches the ~ 3% level. 

It is interesting to consider possible correlations be- 
tween Ag and more usual cosmological parameters. In 
Figures |U and |^ we plot the 1 and 2 cr's confidence level on 
the ns-Xc and w^-Ag planes respectively. As we can see 
there is a strong degeneracy between these parameters. 
Increasing (decreasing) G would yield higher (lower) val- 
ues of ns and lower (higher) values for u>b more consistent 
with CMB data. 

The degeneracy with the scalar spectral index is clear 
since increasing Ag delays recombination, damping the 
small angular scale oscillations. This effect could be 
counterbalanced by increasing ns and the small scale 
power of primordial perturbations. This will also change 
the relative amplitude between odd and even peaks, af- 
fecting the constraints on the baryon density. 

As already described in 0, another possible degener- 
acy is present with the running of the spectral index a s . 
We have therefore considered an extra analysis including 
possible variations in a s . Considering the WMAP data 
only we found Ag = 0.96 ± 0.19. 

Experiment Constraints on Ag at 68% c.l. 



WMAP 


1.01 ±0.16 


WMAP+POL 


0.97 ±0.13 


WMAP+ ACBAR 


1.03 ±0.11 


WMAP+BBN 


0.98 ±0.03 


PLANCK 


1.01 ±0.015 


CVL 


1.002 ±0.004 



TABLE I: Constraints on Ag from current WMAP and 
BBN observations and future constraints achievable from the 
Planck satellite mission and from a cosmic variance limited 
experiment. 





FIG. 4: 68% and 95% likelihood contour plots on the A G - 
ns plane using present CMB data with and without BBN 
constraints. 
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FIG. 5: 68% and 95% likelihood contour plots on the \g- 
uib plane using present CMB data with and without BBN 
constraints. 

As we can sec from Tabic |. future experiments can 
substantially improve the current constraints on Xq. The 
PLANCK Surveyor mission is expected to provide con- 
straints at the ~ 1.5% level. As already discussed in 0, 
the inclusion of polarization data is crucial in breaking 
the degeneracy between Xq and inflationary parameters 
as n s and a s ; we found that neglecting polarization data 
from Planck yields weaker constrains by a factor of ~ 4. 
Polarization data are therefore extremely useful in con- 
straining Xq. 

The ultimate constraint achievable by a cosmic vari- 
ance limited experiment is 0.4%, competitive with cur- 
rent laboratory bounds. 



B. Varying G with redshift 

Here we consider possible constraints on G allowing for 
variations in redshift. Using the simple parameterization 
in Equation U we found that the WMAP data alone yields 
the constraint Xq = 1.01±0.1 at 68% c.L. This constraint 
is better by ~ 40% respect the corresponding bound ob- 
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FIG. 6: 68% and 95% likelihood contour plots on the Xg-zt 
plane using present CMB data. 

taincd with constant G. The reason is due to the extra 
ISW effect that increases the large angular scale CMB 
spectra, in disagreement with the WMAP observations. 

We have then considered a redshift dependence as in 
Equation || with a flat prior < zt < 2. In Figure ^ 
we plot the 68% and 95% confidence levels on the Xq-zt 
plane using only the WMAP data. As we can see, for 
larger values of zt the constraints on Xq are stronger. 
Again, the presence of the ISW effect, irrelevant for zt ~ 
but sizable for larger values, helps in constraining Xq- 

V. CONCLUSIONS 

In this paper we have updated the constraints from 
current CMB data on Newton's gravitational constant 
G. We have found no evidence for variation in this con- 
stant with a constraint of Xq = 1.03 ± 0.11 at 68% c.L 
from WMAP+ACBAR (A G = 0.98 ± 0.03 when BBN 
data is considered). BBN plays therefore a crucial role 
in constraining G. However, even without considering 
the possibility of systematics in current observations of 
primordial elements, the BBN constraints relies on the 
perfect knowledge of the amount of relativistic degrees 
of freedom . Since g* = 5.5 + \N*H any possible extra 
background of relativistic particles, parameterized by the 
effective number of neutrino species N^ff would drasti- 
cally change the BBN bound. Moreover, CMB and BBN 
probe completely different physics and epochs. While 
the agreement between the two results is reassuring, it 
is clear that it would be preferable to have an improved 
and independent CMB constraint. 

We have then considered the constraints achievable 
from ongoing and future satellite experiments. For the 
Planck Surveyor satellite mission we have found a future 
constraints of the order of 1.5% using only CMB data. 
Next, cosmic variance limited experiments as, for exam- 
ple, the future EPIC satellite proposal (see p4|), could 
probe Newton's constant with a ~ 0.4% precision, i.e. 
with grossly the same accuracy currently reached from 
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local experiments. 

It is important to stress that the accuracy on 
achievable by the CMB is limited by how precisely we 
treat the recombination process. Current recombination 
codes should be accurate enough for the Planck mission 
(see e.g. [p5[ ) but this may provide an intrinsic limit 
for the next, beyond Planck, CMB surveys. Moreover, 
recombination could be modified by non-standard mech- 
anisms as dark matter decay or variations in the fine 
structure constant a. High frequency measurements of 
the black-body CMB spectrum, where recombination ab- 
sorption lines are expected, could be helpful in disentan- 
gling the two effects. However, galactic foregrounds at 
those frequencies largely dominate over the CMB signal. 

In this paper we followed a conservative approach by 
considering only future CMB data. It is clear that the in- 
clusion of complementary cosmological data, as expected 
from future galaxy, weak lensing and 21cm surveys, will 
further break the degeneracies between the parameters 
and substantially improve the constraints. We plan to 
discuss this in more detail in a future paper ( ]36| ) . 

Finally, we have considered a variation of G with red- 
shift, parameterizing its variation either with a smooth 
transition between G and XqG, or with a simple step 



function at a transition rcdshift zt- The ISW effect 
arising from redshift variations in G is at odds with the 
low CMB quadrupole measured by WMAP and therefore 
yields stronger constraints on Xq. Current data, also in 
this case, do not exhibit a deviation from the standard 
value. In this respect, the constraints obtained under the 
assumption of constant with redshift G could be consid- 
ered as more conservative. 

Future cosmological data will therefore substantially 
improve the bounds on G and on its possible variations 
with time, space and redshift. By comparing local and 
cosmic measurements, the Newton's constant will be less 
elusive and may shed light on the late accelerated evolu- 
tion of the universe. 
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